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Abstract. In this paper we consider feedback control algorithms for the rapid 
purification of a bipartite state consisting of two qubits, when the observer has access 
to only one of the qubits. We show 1) that the algorithm that maximizes the average 
purification rate is not the same as that that for a single qubit, and 2) that it is always 
possible to construct an algorithm that generates a deterministic rate of purification 
for both qubits. We also reveal a key difference between projective and continuous 
measurements with regard to state-purification. 
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Although quantum measurements are often treated as instantaneous, in practice 
measurement timescales can be significant when compared to other relevant timescales. 
Recent experimental advances have meant that it is now possible to observe continuous 
measurements of individual quantum systems. The measurement of an observable is not 
instantaneous, but takes place over a period of time [Tj [2, EJ E] . This type of continuous 
measurement can be modeled by considering a series of projective measurements on an 
auxiliary system that is weakly coupled to the system of interest. The auxiliary system, 
an environmental degree of freedom, is then averaged out. This produces a continuous 
measurement record which contains information about the evolution of the quantum 
system of interest. The measurement record is then used to construct a best estimate of 
the underlying evolution - which is referred to as an 'unraveling' of the master equation 
for the system EJ El El E] • 

When continuous weak measurement is applied, it is possible to modify the 
evolution via Hamiltonian feedback, where the Hamiltonian applied to the system 
depends on the measurement record [3 El El EH EED, 121 E31 dU [15] . Hamiltonian feedback 
during measurement not only affects the final state of the system, but it can also also 
affect the rate of state reduction. In a protocol described by Jacobs, the average rate 
of state reduction (as measured by the purity of an initially mixed state) for a single 
qubit is maximized [16J. This process is known as rapid state reduction [15J, or as rapid 
purification [163 . Jacobs' protocol is deterministic, but other protocols exist which are 
stochastic and minimize the average time for a single qubit to reach a given purity [17J. 

In this paper we consider the analogous situation of performing rapid state reduction 
in a two qubit system (shown in Fig. [1]). There are two parties, identified as Alice and 
Bob, who may be separated and are not required to communicate. One observer (Alice) 
has access to one qubit, which she can measure and manipulate using local Hamiltonian 
feedback. She does not have access to the second qubit, which is controlled by Bob. This 
corresponds to a physically realistic situation in which Bob's qubit is either spatially 
separated from Alice, or - for architectural reasons - it is not possible to measure Bob's 
qubit directly. Two qubit systems are important because they form the basis for many 
current applications in quantum information processing, and are the simplest system 
which exhibits entanglement. 

In this paper we highlight three aspects of rapid state reduction for bipartite qubit 
states: (1) We highlight a key difference between projective measurements and weak 
measurement. (2) We show that the measurement which provides the maximum rate 
of state reduction for the unprobed qubit in a two qubit system is not necessarily the 
same as that for either of the known optimal single qubit protocols [TBI ITT] . (3) We 
show that it is always possible to purify both qubits deterministically at the same time. 

The bipartite system which Alice and Bob share is described by the density matrix 
p. Without loss of generality, Alice's system undergoes a constant weak measurement, 
giving a measurement record, r(t). It is natural to expand p in the two-qubit Pauli 
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Figure 1. Conceptual layout of the qubits. Alice and Bob both have a single qubit. 
The bipartite system they share is described by the density matrix p. Alice can make 
weak measurements on her system, but not Bob's. 



basis, 

P = 7 E ® a V (!) 

where are real (since p is Hermitian) and lie between —1 and 1, Oj are each of the 
three Pauli matrices and the identity. Each can be found, = Tr pa. t ® <jj. 

The stochastic master equation (SME) which governs the evolution of the density 
operator p in the presence of a weak measurement of a Hermitian observable, y, is given 
by 



dp = -k[y, [y, p]]dt + V2k(yp + py - 2(y)p)dW. (2) 

Here k is the measurement strength. The first term in this equation describes the familiar 
drift towards the measurement axis. The second term in the equation is weighted by 
dW, a Wiener increment with dW 2 = dt. This term describes the update of knowledge 
of the density matrix conditioned on the measurement record |14j . 

Throughout this paper we will consider measurements on Alice's qubit alone. A 
measurement of Alice's qubit along a given axis n means that y is given by 

y = rl ■ al = n x XI + n y YI + n z ZI. (3) 

Here, X, Y and Z are the Pauli matrices, and I is the identity. The tensor product is 
implied. The measurement direction, n may change at each timestep. This is equivalent 
to the application of single qubit Hamiltonian feedback to Alice's qubit, except that the 
measurement axis rotates, rather than Alice's Bloch vector. 

Consider the evolution according to the SME as Alice's qubit undergoes weak 
measurement along the fl axis. The SME (Eqn. ([2])) can be expressed in terms of the 
real Pauli coefficients of the density matrix. If h, m\ and rri2 are mutually orthogonal 
unit vectors, then the corresponding stochastic master equation is given by 

dr mj = - (Akdt + r nI \/8kdW)r mj , (4a) 
dr nj = (r 7i - r nI r nj )V8kdW, (4b) 
dr!j = (r nj - r nI r I:j )VskdW, (4c) 
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for m = {rhx, 777.2} and j = {X, Y, Z, I}. Here, r n j is given by 

r n j = n x r Xj + n y r Yj + n z r Zj , (5) 

and similarly for rh. This is a system of 16 stochastic differential equations. 

In this paper, we will be particularly concerned with the evolution of both Alice 
and Bob's reduced density matrices. Alice's reduced density matrix is given by 

rp ( s I + r XI X + r YI Y + r ZI Z 
p A = Tt b (p) = . (6) 

Bob's reduced density matrix can be described by a similar equation with the indices 
swapped. 

The purity of a quantum state, p can be quantified by the purity, P(p) = Tip 2 . 
Purity has a minimum value of 1/d where d is the dimension of p, and a maximum value 
of 1. The purity of Alice's reduced density matrix is given by 

Pa = TtpI = ^1 + t\), (7) 

where is the Bloch vector of Alice's reduced density matrix. A similar expression can 
be given for Bob's reduced density matrix. 

Consider the rate of state reduction of Alice's qubit as she measures on her own 
system. The change in average purity of Alice's reduced density matrix given in Eqn. 
([7]) under the evolution of the SME, is given by 

dP A = (1 - r 2 nI )(l - r\)Akdt + r^l- r 2 A )V8kdW. (8) 

This expression does not depend on the state of Bob's qubit, and not surprisingly, it is 
identical to the one qubit case. Hamiltonian feedback can be used to implement either 
of the known single qubit Hamiltonian protocols without modification [THJ [T7] . 

Now we consider the opposite situation, when Alice would like to find out the state 
of Bob's qubit. Bob's qubit is not being directly measured. It is only through the 
correlations in the initial system, p, that Alice can learn the state of Bob's system. For 
the most effective purification, Alice and Bob's share an known initial entangled state. 

The average change in purity (according to Alice) of Bob's qubit when measured 
along the ri-axis is given by 

dP B = [(r nX - r nI r IX ) 2 + (r nY - r nI r IY ) 2 + (r nZ - r nI r IZ ) 2 ]AMt + 

[rix(r nX - r nI r IX ) + r IY [r nY - r nI r IY ) + r IZ (r nZ - r nI r IZ )}V8kd^) 

which can be written much more easily by identifying 

Tni = (r n x,r nY ,r nZ ), (10a) 
AR n = r nj - r n/ r 7i . (10b) 

AR n is a vector with three components, one for each of j = {X, Y, Z}. Then the change 
in average purity of Bob's system 

dP B = AR 2 n Akdt + (r E - AR n )V8kdW. (11) 
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We can write the nonlocal correlations matrix: 

Cji = r {j - TuTij, (12) 

for each i,j = {X,Y,Z}. To retrieve AR n from C we multiply by the corresponding 
unit vector, n. That is, 

ARn = Ch. (13) 

The magnitude square of R n gives the rate of state reduction of Bob's qubit, as seen by 
Alice, when she makes a weak measurement along the n direction. 

We now consider how to maximize this rate, in a direct analogy to the optimizing 
strategy used by Jacobs for a single qubit [16J. In particular we find the largest average 
increase in purity of Bob's reduced density matrix. Whilst there are some circumstances 
in which this strategy will not give a globally optimal solution, numerical simulations 
verify its use in this application. The average increase in purity of Bob's reduced density 
matrix is proportional to | Ai? n | 2 . Expressing this in terms of the matrix, C, we wish to 
find 

dP B max) =4kmax\Ch\ 2 . (14) 

h 

The maximum value is given by the largest singular value, a 2 , of C, therefore giving 
a maximum rate increase in purity of Aka\. The value of n which corresponds to this 
maximum rate of state reduction is given by the first right singular vector, v±, of C. 

Point (1) of this paper contrasts weak measurement and projective measurement, 
using a specific example. The measurement which gives the greatest increase in purity 
of Bob's qubit for a projective measurement is not the same as the measurement which 
gives the greatest rate of increase in purity of Bob's qubit for a weak measurement. 

Consider the state, 

(15) 

where — 1 < /3 < 1. is the coherent superposition, with both qubits aligned in the 
z-basis. If P = this is a maximally entangled state, and is entangled unless j3 = ±1. If 
the system of interest undergoes a dephasing process, as is common in many quantum 
systems, then the off diagonal terms of the density matrix decay. If the amount of 
dephasing is characterized by S then the density matrix becomes 

Pt = \ {II + ZZ + (3ZI + pIZ + 7 (IX - YY)) , (16) 

where 7 = a/1 — P 2 — 8- 

The maximum purity from a projective measurement is when Alice measures along 
the z-axis. Any projective measurement by Alice along the z-axis will project Bob's 
state into a pure state (either |0) or |1)), with purity, Pg = 1. As Alice rotates the 
measurement axis, the purity of Bob's state reduces. The matrix C is 

7 

-7 . (17) 

1 - p 2 




C 
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Figure 2. Purification rate of pr, with j3 = 0.5, due to weak measurement of Alice's 
qubit for every orientation of the measurement axis after a dephasing of 6 = 0.01. 4> is 
the zenith angle, and 9 is the azimuthal angle. 



This matrix is diagonal. For small values of dephasing, 5, the maximum rate of state 
reduction occurs in the xy-plane, as is shown in Fig. [2j Therefore the projective 
measurement which gives greatest purity is different from the weak measurement which 
gives the greatest rate of increase in purity. 

Point (2) of this paper contrasts weak measurement for bipartite systems, and 
existing single qubit protocols. Naively, one might expect that the best way for Alice to 
purify Bob's system is to purify her own system fastest by applying known single qubit 
protocols to her own qubit. We show, by giving a specific counter-example, that this is 
not the case. 

Consider the state, 



Pt 




—=(XI + XZ + ZZ) 
v/5 



(18) 



This state is perfectly aligned for Jacobs' one qubit feedback protocol; Alice's reduced 
density matrix lies in the xy-plane [16]. However, for the state pt, the matrix C is given 
by 

(19) 



C 










* 
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Figure 3. Graph showing the rate of average state reduction of Bob's qubit for all 
possible orientations of the measurement axis, and a state given by p t ■ <fi is t ne zenith 
angle, and 9 is the azimuthal angle. 



which has a single non-zero singular value of a% = 1/V5, and a corresponding vector 
of ni = 77| (£ + z). We therefore expect the fastest rate of state reduction for Bob's 
qubit occurs when Alice measures her own system at 45° to the z-axis and to the x-axis, 
as shown in Fig. [3j The numerical values for this plot were k = 0.1, dt = 0.1 and 
iV = 10, 000 repetitions. This measurement does not correspond to either the Jacobs' 
single qubit scheme [IE], or the scheme proposed by Wiseman and Ralph [17]. Purifying 
the unprobed qubit is dependant on the correlations and the nature of those correlations. 
As we show here, it is purified fastest by choosing a measurement axis to make use of 
those correlations, and not by simply using a one qubit protocol. 

Point (3) of this paper is that it is possible to make bipartite state reduction 
deterministic. Increasing the purity of systems deterministically is a desirable property. 
From a theoretical standpoint it makes the equations easier, but the main advantages 
are practical. If the evolution of purity is deterministic, then each qubit is guaranteed 
to reach a set target purity in a given time. When weak measurement is being used 
for state preparation in a multiple qubit system, they will all reach the target purity 
together. 

The condition required for the purity of Alice's reduced density matrix to evolve 
deterministically is that 

r A ■ h = 0. (20) 

That is, the measurement should be taken in the plane orthogonal to the direction 
described by Alice's reduced density matrix. 

The condition for Bob's system to be deterministic is clear from Eqn. (TTTT) . It is 
given by 

r' B -n = 0, (21) 
where r' B = C T r B . Therefore for the evolution of the purity of Bob's reduced density 
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matrix to be deterministic, the measurement axis must be chosen in a plane orthogonal 
to r' B . This plane is spanned by pi and p 2 . 

Deterministic state reduction of Bob's qubit is desirable, but we would still like to 
find the maximum rate of state reduction. Although the measurement axis, n can be 
chosen anywhere in the plane, not all orientations of n will purify Bob's reduced density 
matrix at the same rate. The maximum deterministic rate of state reduction of Bob's 
qubit is found by taking the singular value decomposition of CP P , where P p is given by 



The fastest rate of average deterministic purification of Bob's reduced density 
matrix is given by Akaf where o~\ is the largest singular value of the product CP P , 
and the axis of measurement is given by fi^et = v±, the first right singular vector of CP P . 

It is possible to choose a weak measurement which Alice can make on her qubit 
so that the purity of both her reduced density matrix, and also Bob's reduced density 
matrix both evolve deterministically. If the measurement axis is chosen to be 



then the evolution of the purity of both Alice and Bob's reduced density matrices is 
deterministic. 

In this paper we investigated the effect of weak measurements on a bipartite system 
consisting of two qubits - Alice's qubit and Bob's qubit. We allowed measurement on 
Alice's qubit alone. We gave expressions for the rate of state reduction of both systems, 
based on the measurement record. We have shown how to maximize the average rate of 
state reduction of either system, and how to achieve deterministic evolution of the purity. 
We have demonstrated that weak measurement of two qubits can be very different from 
both projective measurements and the weak measurement of a single qubit. Many 
interesting effects occur in bipartite quantum systems under measurement. 
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